Realization of a Resonant Fermi Gas with a Large Effective Range 
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We have measured the interaction energy and three-body recombination rate for a two-component 
Fermi gas near a narrow Feshbach resonance and found both to be strongly energy dependent. Even 
for deBroglie wavelengths greatly exceeding the van der Waals length scale, the behavior of the 
interaction energy as a function of temperature cannot be described by atoms interacting via a 
contact potential. Rather, energy-dependent corrections beyond the scattering length approximation 
are required, indicating a resonance with an anomalously large effective range. For fields where the 
molecular state is above threshold, the rate of three-body recombination is enhanced by a sharp, 
two-body resonance arising from the closed-channel molecular state which can be magnetically tuned 
through the continuum. This narrow resonance can be used to study strongly correlated Fermi gases 
that simultaneously have a sizeable effective range and a large scattering length. 

PACS numbers: 34.50. Cx, 37.10.Pq, 51.10.+y, 67.85.Lm 



Experimental studies of unitary Fermi gases have ex- 
clusively used energetically broad Feshbach resonances to 
realize strong interactions flT-fli!]. For broad resonances, 
which exhibit negligible energy dependence over the scale 
of the Fermi energy Ep , pairwise interactions in a two- 
component Fermi gas can be completely parameterized 

7H13|. Thus, the equa- 
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by the scattering length, a 
tion of state (EOS) of the Fermi gas in the unitary limit 
(\a\ — > oo) must be a universal function (dependent only 
on Ep for each state and the reduced temperature) that 
describes any dilute Fermi gas with resonant, zero-range 
(i.e. energy- independent) interactions fbH fTa. I20L EH] . 

More generally, however, strongly-interacting Fermi 
systems can also result from resonant, short-range in- 
teractions that are energy dependent and this energy 
dependence can significantly alter the collective behav- 
ior [22-24 [. At minimum, the EOS in the limit \a\ — > 00 
now depends on an additional parameter, the effective 
range r e ff, which characterizes the first-order term in the 
low-energy expansion of the s-wave phase shift 5 (i.e. 
k cotS = —1/a + 2 r eS k 2 + ...). This correction is per- 
tinent to neutron matter at densities in neutron stars 
for which r e g kp <~ 1 where kp is the Fermi wavenum- 
ber 



estimates suggest that the unitary Fermi gas becomes 
even more strongly interacting as |r e ff| is increased from 



zero 
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26] , potentially yielding superfluids with the 
highest critical temperature ever achieved. Most dramat- 
ically, the very existence of certain novel phases of matter 
(e.g. the breached-pair superfluid phase [27J) depend on 
interactions having a particular momentum dependence. 

Narrow Feshbach resonances (FR) should allow for the 
study of such behavior in an atomic Fermi gas as they 
have been predicted to exhibit a strong ener gy d epen- 
dence over the scale of E F [T3-ES 0, 0, EsQ- The 
effective range r c g in the vicinity of a FR is inversely pro- 
portional to its width and, for narrow resonances, can be 
large in comparison to 1/kp or even the van der Waals 



length 4dw [19| . potentially allowing for the study of 
dilute Fermi gases with large a and large r e g. 

In this Letter, we experimentally demonstrate that in- 
teractions in a two-component Fermi gas near a narrow 
FR cannot simply be parameterized by a field-dependent 
scattering length but require an energy-dependent scat- 
tering length or, near the resonance, an anomalously 
large effective range. This is done by comparing the mea- 
sured interaction energy for a two-component Fermi gas 
at different temperatures with predictions of mean-field 
theory. The interaction energy is measure d by radio- 
frequency spectroscopy using a third state KBfJ|. Mea- 
surements with cold gases accurately determine the lo- 
cation and width of the resonance. For hotter samples, 
the field-dependence of the interaction energy exhibits an 
asymmetric lineshape which we interpret to be a signa- 
ture of energy dependent interactions. We further show 
that inelastic loss from the gas by three-body recombi- 
nation is enhanced by a sharp resonance in the two-body 
continuum as the molecular state is tuned above thresh- 
old. From the inelastic loss data we can extract the 
vibrational relaxation rate constant due to atom-dimer 
collisions for the molecular state associated with the res- 



25 . In other cases where r e g < 0, theoretical onance. 



In our experiments, we investigate the narrow reso- 
nance which occurs in s-wave collisions between the two 
lowest-energy hyperfine states in 6 Li (states |1) and |2)) 
near a field of Boa = 543.3 G [3l[. Long-lived 6 Li2 
molecules had previously been created using this reso- 
nance 32j. In the vicinity of the resonance, the scatter- 
ing length varies as a(B) = ab g (1 — A/(£? — Boo)) (see 
Fig. HT.a)). Here, A = 0.1 G is the resonance width 
and dbg — 62 clq is the background scattering length 
which varies relatively slowly with field. This narrow 
FR occurs when the least bound vibrational level of 
the X 1 E+ singlet potential crosses threshold. The en- 
ergy of the molecular state relative to threshold varies 
as E c = \i r (B — Boo) where [i r is the difference between 
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FIG. 1. Narrow Feshbach resonance in 6 Li. (a) The scattering 
length as a function of magnetic field, (b) For E < 0, the 
binding energy of the molecular state relative to threshold 
is shown. The molecular state crosses threshold at Boo = 
543.286(3) G. For E > 0, a contour plot of sin 2 S(E, B) 
demonstrates the sharp resonance that the molecular state 
creates as it is tuned through the continuum, (c) Fraction 
of atoms remaining after a 100 ms hold time at various B- 
fields when the resonance is approached from below (o) or 
above (•) showing an asymmetric loss resonance for a cloud 
at T = 7fJ,K. 



the magnetic moment of the colliding atoms and the bare 
molecular state given in this case by \i r = 2 \i& where p,B 
is the Bohr magneton. 

Near threshold, the s-wave scattering phase shift S is 



5bg 



tan 



r/2 

E-E c 



(1) 



the sum of a background phase shift o\, g = — k Ob g and 
a resonant phase shift with a Breit-Wigner form [l9| . 
Resonance occurs when the relative kinetic energy E — 
h 2 k 2 /m is equal to the energy of the molecular state rela- 
tive to threshold E c . Here, hk is the relative momentum 
and m is the mass of a 6 Li atom. Near threshold, the 
energy dependent resonance width T(E) = 2 fc Ob g /i r A 
describes the finite lifetime of the bound state due to 
coupling to the continuum. A contour plot of sin 2 5 in 
the upper half of Fig.QJb) illustrates the sharp resonance 
that tunes through the continuum when the molecular 
state is above threshold. 

The phase shift in Eq. [T] determines the scattering am- 
plitude / = l/(k cot 5 — ik) where, as in Ref. [26|, we 



h A k A 



k cot S = 

dbg 



\i r (B — Boo ) 



K 2 k 2 



fir {B - Boo) + yUr A ' 
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assuming that k a^ g -C 1 ■ Note that near the resonance 
(i.e. for (B — B^) <C A) and near threshold, the scat- 
tering phase shift is given by the effective range expan- 



sion, k cot 5 = — — 



and 



r e g = 



2h 2 



cff k 2 where a ros = a bg 



For the narrow FR studied 



n a hs fi r A ' 

here, the effective range, being inversely proportional to 
(dbg A), has an anomalously large value r e g = — 7x 10 ao 
with a magnitude much larger than the average inter- 
particle spacing in the dilute Fermi gases studied here. 
Narrow FRs in 6 Li- 40 K mixtures should also exhibit 
l r off| 3* f v dw, though the associated energy dependence 
has not been verified 

To study the narrow FR in 6 Li we have (1) monitored 
the loss of atoms as a function of magnetic field, and 
(2) performed radio-frequency (rf) spectroscopy to mea- 
sure interaction energy shifts using the transition be- 
tween state |2) and the third lowest energy hypcrfine 
state in 6 Li (state |3)). We can prepare degenerate or 
thermal mixtures of 6 Li atoms with equal populations 
in states |1) and |2) b y e vaporative cooling in a single 
or crossed-dipole trap j35||. To control the density and 
temperature of atoms in our experiments, we adjust the 
degree of evaporative cooling and the final trap depth 
following evaporation. 

Our initial observation of the narrow FR is shown in 
Fig. [ljc) which plots the fraction of atoms remaining 
in the optical dipole trap after a 100 ms hold time at 
each field. For this data, the average density per spin 
state n ~ 10 13 cm -3 , the temperature T ~ 7/iK, and 
the o's (•'s) show the fraction remaining when the FR 
is approached from a field well-below (well-above) reso- 
nance at a ramp-rate ~ 50mG/ms. Atom loss from this 
system must be due to three-body recombination since 
two-body inelastic processes are frozen out at this tem- 
perature for a mixture of 6 Li atoms in states |1) and 
1 2). This feature is asymmetric with loss preferentially 
occurring on the high-field (atomic) side, in contrast to 
the behavior observed for Fermi gases near broad reso- 
nances (36J. We show that the asymmetry results from 
the fact that the molecular state enhances three-body re- 
combination as it tunes through the continuum. Higher- 
temperature gases can have collisions sufficiently ener- 
getic to satisfy the resonance condition E = fi r (B — Boo) 
even for magnetic fields far-detuned above resonance (i.e. 
for k B T - fi r (B - Boo) > Mr A). Thus, the width of 
the loss feature on the high-field side of resonance grows 
as the temperature is increased. 

To investigate the dependence of loss on temperature 
and field near this FR, we monitor the number of atoms 
as a function of time at fixed fields with samples initially 
prepared above resonance at one of two different temper- 
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FIG. 2. Recombination loss rates for gases initially prepared 
at (a) Ti = 3.2 /iK and (b) T< = 10 /iK. The solid lines are 
single-parameter fits to a model (Eq. [4]) which assumes that 
the molecular state resonantly enhances loss as it is tuned 
through the continuum. 



atures. The density decays according to h = —L3 n 3 due 
to three-body recombination where L3 is the three-body 
recombination loss rate constant. For initial Tj = 3.2 /iK 
(Ti = 10 /iK), we have extracted L3 for different B-ficlds 
as shown in Fig. [Ha) (Fig. [5] (b)). In both cases, the 
sample is initially at T ~ 0.5 Tp, allowing us to analyze 
this weakly degenerate system as a thermal gas [35| . As 
the B-field approaches resonance, the temperature of the 
cloud increases as shown in the Fig. [2] insets. A smooth- 
ing spline is fit to the temperature data (solid lines in 



Fig. [2] insets) and used to compute n and extract L 3 [37 1. 

To understand the behavior of L3 above resonance, we 
consider the fact that two atoms with relative kinetic en- 
ergy E ~ E c are resonantly coupled to the vibrationally 
excited molecular state responsible for the FR which can 
decay in a subsequent atom-dimcr collision. A narrow 
FR which couples two colliding atoms to a resonant state 
which has a finite lifetime (H/Tq) due to inelastic decay, 
induces loss with a two-body loss rate coefficient [38| 



x 2 = AI 



T(E)T 



m k (E-2fi B (B - Boo))' 
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Here, both inelastic loss and coupling to the continuum 
contribute to determining the total linewidth r to t = To + 
T(E) . We assume that r = h K ad (ni + n 2 ) where if at j 
is the atom-dimer relaxation rate coefficient and n\ (112) 
is the density of atoms in state |1) (|2)). For each state, 
the density decays according to h = — L^n 3 since L2 is 
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FIG. 3. Rabi spectra of the |3) — > |2) transition for a 95 fjs- 
duration 7r-pulse. (a) With no atoms in state |1), the spec- 
trum accurately determines the B-field. (■) With atoms in 
state |1) present, the lineshape is broadened and shifted by 
interactions. The data is well-approximated (dashed) by con- 
volving the Rabi transition lineshape (dot-dashed) with the 
lineshape obtained for a thermal gas in a harmonic trap as- 
suming a frequency shift proportional to the local density. 
Fits which allow the Rabi and offset frequencies to vary (solid 
line) are used to determine the mean shift v. 



itself proportional to density. A simple expression for 
L3 is obtained after performing a thermal average if we 
assume that the loss resonance is sharply peaked (i.e. 
Ttot <C ksT) and the field is detuned from resonance 
such that T(fi r (B — Boo)) ^ IV Under these conditions 



U = 



2h 3 K ad 
(7rmfc B T) 3 / 2 



exp 



/i r (B — Boo) 



knT 



(4) 



which has a 1/e width given by ks T//j, r . 

We have fit both the 3.2 /iK and 10 /iK data sets using 
Eq. [4] with K ad as the only free parameter (solid lines in 
Fig. [2]). To evaluate L 3 in Eq.gJ we use Boo = 543.286 G 
(determined from interaction energy shifts described be- 
low), the experimentally measured temperature (insets 
in Fig. [5]), and the unknown atom-dimer relaxation rate 
constant K^. Eqn. [4] provides an excellent descrip- 
tion of the 3.2 ijK data and a good description of the 
10 /iK data for atom-dimer relaxation rate coefficients 
A' ad = 3(1) x 10~ 10 cm 3 /s and A ad = 2(1) x 10- 10 cm 3 /s 
respectively. Here, the uncertainty in A a( j is due to 
the uncertainty in trap frequencies. From the aver- 
age of both fits and including our systematic uncer- 
tainty in atom number (~ 30%) we determine A a( j = 
2.5 ± 1.9 x 10 _10 cm 3 /s. This narrow resonance can po- 
tentially be used as an energy-selective "knife" for forced 
evaporative cooling j38|. A slight reduction in tempera- 
ture observed in both data sets when \x r (B-B m ) > 2 Ep 
(see Fig. [5] insets) may be a result of such evaporation. 

To measure the dependence of interaction shifts on 
field and temperature, we have performed rf spectroscopy 
between states |3) and |2) (Fig. [3]). The presence of an 
incoherent population in state |1) shifts the |3) — > |2) 
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transition frequency by the difference in interaction en- 
ergy for a state- 1 2) or a state-) 3) atom interacting with 
the background density of atoms in state |1) [30|. The 
two-body contribution to mean-field interactions is de- 
termined by the real part of the complex scattering am- 
plitude, /, yielding a shift of the rf transition given by 



Ai/ 



2h 



(-Rc[/ 13 ]+Rc[/ 12 ]) ni 



(5) 



where ri\ is the density of atoms in state |1). Inter- 
actions between atoms in states |1) and |3) can be de- 
scribed by a contact potential parameterized by a scat- 
tering length ai3 for which — Re[/i3] = ai3 at low energy 
(i.e. kai3 <C 1). However, the narrow |1) — | 2} - resonance 
produces a scattering amplitude with a real part that is 
energy dependent. In this case, /12 = l/(fc cot £12 — ik) 
where k cot 5\2 is given by Eq. [5] For a thermal gas at 
temperature T with density ri\, the expected frequency 
shift is given by Eq. [5] where Re[/i2] is replaced by its 
thermal average to yield a T-dependcnt frequency shift 
Av(rii,T, B). Since the density distribution is inhomo- 
geneous in the trap, the rf-transition lineshape is both 
shifted and broadened [1^ relative to that observed when 
atoms in state |1) are absent (Fig. [3]). In order to avoid 
inelastic loss and heating while adjusting B, we prepare 
a 50/50 mixture of atoms in states |1) and |3) (rather 
than |1) and |2)), then shift to the field of interest and 
finally drive the rf transition from |3) — > |2) immediately 
before absorption imaging. 

In Fig. |4] we plot — 3% = for gases at different absolute 
temperatures where v is the mean value of the observed 
lineshape (relative to the bare transition frequency for 
an isolated atom) and n is the average density. In terms 
of the mean-field prediction, the plotted quantity equals 
a* = (-Rc[/i 2 (a hg ,B,B 00 , A)]) T - ai 3 . Here (...) T in- 
dicates a thermal average. The lowest temperature data 
(Fig. H£a)) is fit to A x a*(T, a i3 , a bg , B, Boo, A) where 
the amplitude A accommodates for our uncertainty in 
the absolute number and trap frequencies. To minimize 
X 2 , the amplitude A, resonance location Boo, and width 
A are varied while the measured temperature T and the 
slowly varying scattering lengths (113(B) ~ — 267 ao and 
a,bg(B) ~ 62 ao are fixed. We find the minimum x 2 
for Boo = 543.286(3) G and A = 0.10(1) G where the 
quoted errors include the uncertainty in field calibration. 

Using Boo = 543.286 G and A = 0.10 G determined 
from the lowest temperature data, we plot A x a* in 
Figs. |U[a-d) for each specified temperature where only A 
is allowed to vary (solid lines). For comparison, the pre- 
diction for a contact potential with — Re[/i2] = ab 9 (l — 
A/(B — Boo)) is also shown (dashed lines). Clearly, using 
a contact interaction to model |1) — |2)-scattering fails 
to describe the observed interaction shifts as the tem- 
perature is increased. The mean-field predictions which 
assume an energy-dependent Rc[/i 2 ] provide a much bet- 
ter description of the data. In particular, the increase of 
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FIG. 4. Interaction-induced frequency shifts for clouds at dif- 
ferent temperatures, (dashed-lines) Shifts predicted if a con- 
tact potential describes |1) — |2) collisions, (solid lines) Shifts 
predicted by a mean-field theory that includes the energy de- 
pendence of the scattering phase shift expected near a narrow 



resonance. 



the lineshape asymmetry with temperature, the behav- 
ior of the interaction energy on the high-field side of the 
resonance and the field at which the strongest attractive 
interaction energy is observed are well-described by this 
model. However, even the energy-dependent model does 
not describe the data for fields extremely close to reso- 
nance (e.g. the T = 10 fiK data in Fig. Be)). These dis- 
crepancies may be due to correlations in the many-body 
wavefunction not described by mean-field theory. A re- 
cent calculation of the interaction energy near a narrow 
FR using a virial expansion (applicable for large T/Tf) 
also predicts asymmetric behavior across the resonance 
but does not quantitatively describe our data [26| . 
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In summary, we demonstrated that the inelastic loss of 
fcrmions near a narrow FR is enhanced when the molec- 
ular state is above threshold, contrary to the behavior 
observed near broad resonances. We also demonstrated 
that interatomic interactions near a narrow FR cannot 
be described by a contact potential. A mean-field theory 
which includes the full energy dependence of the scat- 
tering phase shift expected for a narrow FR provides a 
good, though not perfect, description of the observed in- 
teraction shifts. This FR can be used to study strongly- 
correlated Fermi systems with energy-dependent interac- 
tions, potentially allowing for superfluids with the high- 
est critical temperature ever achieved 23, Hf|, the ob- 
servation of a breached-pair supcrfluid phase [27(, and 



the determination of the equation of state for neutron 
star matter at densities comparable to or higher than 
the neutron drip density [25| . 
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